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INTRODUCTION 

It is generally accepted that QCD with two massless quarks undergoes a chiral 
O \ phase transition (see reviews by DeTar^, Ukawa^ and Smilga^ for recent results on this 

^ \ topic). This leads to important observable signatures in the real world with two light 

■ quarks. The order parameter of the chiral phase transition is the chiral condensate 
^ ' (^^)- Below the critical temperature chiral symmetry is spontaneously broken with 

('?/''?/') 7^ 0. One consequence is that parity doublets are absent in the hadronic spectrum. 
For example, the pion mass and the 5 mass are very different. To a good approximation, 
Q \ as dictated by the Goldstone theorem, the pion is massless. According to lattice QCD 

■ simulations, chiral symmetry is restored for T > where ('?/''?/') = 0. For two light 
flavors the critical temperature is expected to be ~ 140 MeV . In the restored phase, 

^ . parity doublets are present, and a massive pion is degenerate with the cr- meson. It is 

well known that the QCD Lagrangian has two chiral symmetries: the f/A(l) symmetry 
and the SU{Nf) x SU{Nf) symmetry. As was in particular pointed out by Shuryak"^ 
not necessarily both symmetries are restored at the same temperature^. This may lead 
to interesting physical consequences. 

According to the Banks-Casher formula^, the chiral condensate is directly related 
to the average spectral density of the Dirac operator near zero virtuality. However, 
the eigenvalues of the Dirac operator fluctuate about their average position over the 
ensemble of gauge field configurations. The main question we wish to address in these 
lectures is to what extent such fluctuations are universal. If that is the case, they do 
not depend on the full QCD dynamics and can be obtained from a much simpler chiral 
Random Matrix Theory (chRMT) with the global symmetries of the QCD partition 
function. 

This conjecture has its origin in the study of spectra of complex systems^. Ac- 
cording to the Bohigas conjecture, spectral correlations of classically chaotic quantum 
systems are given by RMT. A first argument in favor of universality in Dirac spectra 
came from the analysis of the finite volume QCD partition function*^. As has been 
shown by Gasser and Leutwyler^, for box size L in the range 

l/A<L<l/m^, (1) 
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(A is a typical hadronic scale and is the pion mass) the mass dependence of the 
QCD partition function is completely determined by its global symmetries. As a con- 
sequence, fluctuations of Dirac eigenvalues near zero virtuality are constrained by, but 
not determined by, an infinite family of sum rules^ (also called Leutwyler-Smilga sum 
rules). For example, the simplest Leutwyler-Smilga sum rule can be obtained from 
the microscopic spectral density^'^ (the spectral density near zero virtuality on a scale 
of a typical eigenvalue spacing). On the other hand, the infinite family of Leutwyler- 
Smilga sum rules is not sufficient to determine the microscopic spectral density. The 
additional ingredient is universality. A priori there is no reason that fluctuations of 
Dirac eigenvalues are in the same universahty class as chRMT. Whether or not QCD is 
inside this class is a dynamical question that can only be answered by full scale lattice 
QCD simulations. However, the confldence in an affirmative answer to this question 
is greatly enhanced by universality studies within chiral Random Matrix Theory. The 
aim of such studies is to show that spectral fluctuations do not depend of the details of 
the probability distribution. Recently, it has been shown that the microscopic spectral 
density is universal for a wide class of probability distributions^"^' We will 

give an extensive review of these important new results. 

The fluctuations of Dirac eigenvalues near zero virtuality are directly related to the 
approach to the thermodynamic limit of the chiral condensate. In particular, knowledge 
of the microscopic spectral density provides us with a quantitative explanation^^ of flnite 
size corrections to the valence quark mass of dependence of the chiral condensate^^ . 

Because of the ?7a(1) symmetry of the Dirac operator two types of spectral fluc- 
tuations can be distinguished. Spectral fluctuations near zero virtuality and spectral 
fluctuations in the bulk of the spectrum (Actually, there is a third type: spectral fluc- 
tuations near the end-points of lattice QCD Dirac spectra. However, this region of the 
spectrum is completely unphysical, and it will not be considered in these lectures.) 

Recently, it has become possible to obtain all eigenvalues of the lattice QCD 
Dirac operator on reasonably large lattices^^' making a direct veriflcation of the 
above conjecture possible. This is one of the main objectives of these lectures. This is 
easiest for correlations in the bulk of the spectrum. Under the assumption of spectral 
ergodicity^^, eigenvalue correlations can be studied by spectral averaging instead of 
ensemble averaging^^' On the other hand, in order to study the microscopic spectral 
density, a very large number of independent gauge fleld conflgurations is required. First 
lattice results conflrming the universality of the microscopic spectral density have been 
obtained recently^°. 

At this point I wish to stress that there are two different types of applications of 
Random Matrix Theory. In the first type, fluctuations of an observable are related to 
its average. Because of universality it is possible to obtain exact results. In general, 
the average of an observable is not given by Random Matrix Theory. There are many 
examples of this type of universal fluctuations ranging from atomic physics to quantum 
fleld theory (a recent comprehensive review was written by Guhr, Miiller-Groeling 
and Weidenmiiller^^). Most of the examples are related to fluctuations of eigenvalues. 
Typical examples arc nuclear spectra^^, acoustic spectra^^, resonances in resonance 
cavities^'', ^-matrix fluctuations^^' and universal conductance fluctuations^^. In these 
lectures we will discuss correlations in the bulk of Dirac spectra and the microscopic 
spectral density. The second type of application of Random Matrix Theory is as a 
schematic model of disorder. In this way one obtains qualitative results which may be 
helpful in understanding some physical phenomena. There are numerous examples in 
this category. We only mention the Anderson model of localization^^, neural networks^^, 
the Gross- Witten model of QCD^^ and quantum gravity^"^ . In these lectures we will 
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discuss chiral random matrix models at nonzero temperature and chemical potential. In 
particular, we will review recent work by Stephanov^^ on the quenched approximation 
at nonzero chemical potential. 

At nonzero chemical potential the QCD Dirac operator is nonhermitean with eigen- 
values scattered in the complex plane. As was first pointed out by Fyodorov et al.^^, 
this leads to the possibility of a new type of universal behavior. Characteristic fea- 
tures of Dirac spectra will be discussed at the end of this lecture. For a review of 
nonhermitean random matrices, we refer to the talk by Nowak^'' in these proceedings. 

In the first lecture we will review some general properties of Dirac spectra including 
the Banks-Casher formula. Prom the zeros of the partition function we will show that 
there is an intimate relation between chiral symmetry breaking and correlations of 
Dirac eigenvalues. Starting from Leutwyler-Smilga sum-rules the microscopic spectral 
density will be introduced. We will discuss the statistical analysis of quantum spectra. 
It will be argued that spectral correlations of 'complex' systems are given by Random 
Matrix Theory. We will end the first lecture with the introduction of chiral Random 
Matrix Theory. 

In the second lecture we will compare the chiral random matrix model with QCD 
and discuss some of its properties. We will review recent results that show that the 
microscopic spectral density and eigenvalue correlations near zero virtuality are strongly 
universal. Lattice QCD results for the microscopic spectral density and correlations in 
the bulk of the spectrum will be discussed in detail. We will end the second lecture 
with a review of chiral Random Matrix Theory at nonzero chemical potential. New 
features of spectral universality in nonhermitean matrices will be discussed. 



THE DIRAC SPECTRUM 
Introduction 

The Euclidean QCD partition function is given by 

ZqcBim, e)^ j dAdet {^D + m)e-^^^/'^+^^^ (2) 

where jD = ^^{0^ + iA^) is the anti-Hermitean Dirac operator and Sym is the Yang- 
Mills action. The integral over field configurations includes a sum over all topological 
sectors with topological charge u. Each sector is weighted by exp(i^i/). Phenomenolog- 
ically the value of the vacuum ^-angle is consistent with zero. We use the convention 
that the Euclidean gamma matrices are Hermitean with {7^,7j/} = '26^^. The inte- 
gral is over all gauge field configurations, and for definiteness, we assume a lattice 
regularization of the partition function. 

Our main object of interest is the spectrum of the Dirac operator. The eigenvalues 
Afc are defined by 

-fD(j)k = iXkh- (3) 

The spectral density is given by 

p{X)^Y.^{X-X,). (4) 

k 

Correlations of the eigenvalues can be expressed in terms of the two-point correlation 
function 

P2(A,A') = (p(AKA')), (5) 
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where (■ • ■) denotes averaging with respect to the QCD partition function (^). The 
connected two-point correlation function is obtained by subtraction of the product of 
the average spectral densities 

p,(A,A')=P2(A,A')-(p(A))(p(A')). (6) 

Because of the Ua{1) symmetry 

{75,7^} = 0, (7) 

the eigenvalues occur in pairs ±A or are zero. The eigenfunctions are given by (pk 
and 750fc, respectively. If 'y^cpk = ±0/c, then necessarily A^ = 0. This happens for a 
solution of the Dirac operator in the field of an instanton. In a sector with topological 
charge u the Dirac operator has u exact zero modes with the same chirality. In order 
to represent the low energy sector of the Dirac operator for field configurations with 
topological charge u, it is natural to choose a chiral basis with n right-handed states 
and m = n + u left-handed states. Then the Dirac matrix has the block structure 




where T is an x m matrix. For m = 2 and n = 1, one can easily convince oneself 
that the Dirac matrix has exactly one zero eigenvalue. We leave it as an exercise to the 
reader to show that in general the Dirac matrix has \m — n\ zero eigenvalues. 

In terms of the eigenvalues of the Dirac operator, the QCD partition function can 
be rewritten as 

ZqcdK^) = Ee^'^'n^/' / dAl[l[iXl + m})e"'-^'/^ (9) 

/ f k 

where /j, dA denotes the integral over field configurations with topological charge z/, 
and 11/ is the product over Nj flavors with mass m/. The partition function in the 
sector of topological charge v is obtained by Fourier inversion 

Zu{m) = — / dee-''''Zc,cx^{m, 9). (10) 

ZTT JO 

The fluctuations of the eigenvalues of the QCD Dirac operator are induced by the 
fluctuations of the gauge fields. Formally, one can think of integrating out all gauge 
field configurations for fixed values of the Dirac eigenvalues. The transformation of 
integration variables from the fields. A, to the eigenvalues, A^, leads to a nontrivial 
" Jacobian" . Universality in Dirac spectra has its origin in this " Jacobian" . 



The free Dirac spectrum can be obtained immediately from the square of the Dirac 
operator. For a box of volume Li x L2 x L3 x L4 one finds 




where we have used periodic boundary conditions in the spatial directions and anti- 
periodic boundary conditions in the time direction. The spectral density is obtained 
by counting the total number of eigenvalues in a sphere of radius \L/2'k. The result is 

p'^°^(A) ~ V\\ (12) 



4 



1.0 



0.8 
P(A) 
0.6 



0.4 



0.2 



0.0 

0.0 1.0 2.0 3.0 4.0 5.0 

K 

Figure 1. The free Dirac spectral density (dotted curve) and the spectral density of the 
Kogut-Susskind Dirac operator for a gauge field configuration generated with f3 = 2.4 (histogram). 
Both spectral densities are on a 12^ lattice and are normalized to unit area. 

For future reference, we note that in the generic case, when the sides of the hypercube 
are related by an irrational number, the eigenvalues are uncorrelated, i.e. 

P2(A,A') = (p(A))(p(A')). (13) 

Two examples of Dirac spectra are shown in Fig. 1. The dotted curve represents 
the free Kogut-Susskind Dirac spectrum on a 12'* lattice with periodic boundary condi- 
tions in the spatial directions anti-periodic boundary conditions in the time direction. 
For an Ni x N2 x x N4 lattice this spectrum is given by 

^ 2 (s... . s... (^) . (^) . ... (^))" . (14, 

Here, rii = 0, 1, ■ ■ ■ , [N,/2] {i = 1, 2, 3) and ^4 = 0, 1, ■ ■ ■ , [{N^ + l)/2]. The Kogut- 
Susskind Dirac spectrum for an SU{2) gauge field configuration with j3 = 2.4 on the 
same size lattice is shown by the histogram in the same figure (full curve). We clearly 
observe an accumulation of small eigenvalues. 

The Banks-Casher Relation 

The order parameter of the chiral phase transition, (i/jip), is nonzero only below 
the critical temperature. As was shown by Banks and Casher^, (ipip) is directly related 
to the eigenvalue density of the QCD Dirac operator per unit four-volume 

S = 1(^^)1 = lim^^^. (15) 

It is elementary to derive this relation. The chiral condensate follows from the partition 
function (^ (all quark mass are chosen equal), 

- 1 1 2m 

(^^) = - lim TTT^dm log Z{m) = - lim - Tj— — ^) • (16) 
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If we express the sum as an integral over the average spectral density, and take the 
thermodynamic limit before the chiral limit, so that many eigenvalues are less than 
m, we recover (|T3]). The order of the limits in (|T3]) is important. First we take the 
thermodynamic limit, next the chiral limit and, finally, the field theory limit. As can 
be observed from (|TBp the sign of {ipip) changes if m crosses the imaginary axis. 



An important consequence of the Bank-Casher formula ([T5|) is that the eigenvalues 
near zero virtuality are spaced as 

AA = l/p(0) = tt/SI/. (17) 

This should be contrasted with the eigenvalue spectrum of the non-interacting Dirac 
operator. Then one obtains from (|1^) an eigenvalue spacing equal to AA ~ l/V^^^. 
Clearly, the presence of gauge fields leads to a strong modification of the spectrum 
near zero virtuality. Strong interactions result in the coupling of many degrees of 
freedom leading to extended states and correlated eigenvalues. On the other hand, 
for uncorrelated eigenvalues, the eigenvalue distribution factorizes, and for A 7^ 0, we 
have p(A) ~ A^^^^^ in the chiral limit, i.e. no breaking of chiral symmetry. One 
consequence of the interactions is level repulsion of neighboring eigenvalues. Therefore, 
the two smallest eigenvalues of the Dirac operator, ±Amin repel each other, and the 
Dirac spectrum will have a gap at A = with a width of the order of 1/Ty . 



Spectral Correlations and Zeros of the Partition Function 

The study of zeros of the partition function has been a fruitful tool in statistical 
mechanics'^^' In QCD, both zeros in the in the complex fugacity plane and the 
complex mass plane have been studied^^' Since the QCD partition function is a 
polynomial in m it can be factorized as (all quark masses are taken to be equal to m) 

Zq^cJi{m,d) = W{m-mk). (18) 

k 

Because configurations of opposite topological charge occur with the same probability, 
the coefficients of this polynomial are real, and the zeros appear in complex conjugate 
pairs. For an even number of flavors the zeros occur in pairs im^. In a sector with 
topological charge this is also the case for even Nj x v. The chiral condensate is 
given by 

1 11 
S(m) = - lim TT^dra log ^qcd("^, 0) = - lim • (^9) 

For an even number of flavors, S(m) is an odd function of m. In order to have a 
discontinuity at m = 0, the zeros in this region have to coalesce into a cut along the 
imaginary axis in the thermodynamic limit. 

In the hypothetical case that the eigenvalues of the Dirac operator do not fluctuate 
the zeros are located at = ±iXk- In the opposite case, of uncorrelated eigenvalues, 
the eigenvalue distribution factorizes and all zeros are located at iicx, where = (Af,). 
As a result, the chiral condensate does not show a discontinuity across the imaginary 
axis and is equal to zero. 

We hope to convince the reader that the presence of a discontinuity is intimately 
related to correlations of eigenvalues^^. Let us study the effect of pair correlations for 
one flavor in the sector of zero topological charge. The fermion determinant can be 
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written as 

(n(^^ + AD) = E ( ^ ) ^^^^-'^^A? • • • A^. (20) 

There are 

( 21 ) - D" 

ways of selecting I pairs from • • • A^. The average of each pair of different eigenvalues 
is given by 

(A^,A^) ^a' + C2, (21) 
where is the expectation value of A^ and C2 is the connected correlator 

^2 = (A^A^)-(A^)(A^), m^n. (22) 
This results in the partition function 



N 12 



*1 ' X , 

N \ k 



Z{m) = E E^'^"^-'^ \ k][2l\ - (23) 



fc=0 l=Q 



After interchanging the two sums, one can easily show that Z{m) can be expressed as 
a multiple of a Hermite polynomial 



Z{m) = (-f )^/^H^((a^ + m')/sJ-2C,). (24) 

In terms of the zeros of the Hermite polynomials, Zk-, the zeros of the partition function 
are located at 



ml = ZkyJ-2C2 - (25) 



Asymptotically, the zeros of the Hermite polynomials are given by ~ Tik/2yN. In 
order for the zeros to join into a cut in the thermodynamic limit, they have to be spaced 
as ~ This requires that 

C2 ~ (26) 

The density of zeros is then given by 

dk , , 

Nm. (27) 

dm 

We conclude that pair correlations are sufficient to generate a cut of Z{m) in the 
complex m-plane, but the chiral symmetry remains unbroken. Pair correlations alone 
cannot suppress the effect of the fermion determinant. 

Leutwyler-Smilga Sum Rules 

We have shown that pair-correlations are not sufficient to generate a discontinuity 
in the chiral condensate. In this subsection we start from the assumption that chiral 

symmetry is broken spontaneously, and look for consistency conditions this imposes 
on the Dirac spectrum. As has been argued by Gasser and Leutwyler^ and Leutwyler 
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and Smilga^, in the mesoscopic range (0), the mass dependence of the QCD partition 
function is given by (for simphcity, all quark masses have been taken equal) 

Z'^{m,e)r^f rft/e^^^^^^Trf/e^^/-/^ ^28) 

JU£G/H 

The integral is over the Goldstone manifold associated with chiral symmetry break- 
ing from G to H. For three or more colors with fundamental fermions G/H = 
SU{Nf) X SU{Nf)/ SU{Nf). The finite volume partition function in the sector of 
topological charge u follows by Fourier inversion according to (]T0|). The partition func- 
tion for u = is thus given by (|28|) with the integration over SU{Nf) replaced by an 
integral over U{Nf). The case of Nf = 1 is particularly simple. Then only a U{1) 
integration remains, and the partition function is given by^ Z^^Q^m) = Io{mVT,). Its 
zeros are regularly spaced along the imaginary axis in the complex m-plane, and, in 
the thermodynamic limit, they coalesce into a cut. 

The Leutwyler-Smilga sum-rules are obtained by expanding the partition function 
Zy{m) in powers of m before and after averaging over the gauge field configurations and 
equating the coefficients. This corresponds to an expansion in powers of m of both the 
QCD partition function (Q) and the finite volume partition function (|28| ) in the sector 
of topological charge v. As an example, we consider the coefficients of in the sector 
with z/ = 0. This results in the sum-rule 



(29) 



where the prime indicates that the sum is restricted to nonzero positive eigenvalues. 

The next order sum rules are obtained by equating the coefficients of order m^. 
They can be combined into 

We conclude that chiral symmetry breaking leads to correlations of the inverse eigen- 
values. However, if one performs an analysis similar to the one in previous section, it 
can be shown easily that pair correlations given by (|30D do not result in a cut in the 
complex m-plane. Apparently, chiral symmetry breaking requires a subtle interplay of 
all types of correlations. 

For two colors with fundamental fermions or for adjoint fermions the pattern of 
chiral symmetry breaking is different. Sum rules for the inverse eigenvalues can be 
derived along the same lines. The general expression for the simplest sum-rule can be 
summarized as^^' 



(31) 



XI' 4(|z/| + (dim(G/iJ) + l)/Ar^)- 
The Leutwyler-Smilga sum-rules can be expressed as an integral over the average 



spectral density and spectral correlation functions. For the sum rule (|29|) this results 



m 



1 f {p{X))d\ _ 1 

^2^2 7 A2 4iV/ ^ ' 

If we introduce the microscopic variable 

u = XVE, (33) 
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this integral can be rewritten as 



u ,,du 1 



The thermodynamic hmit of the combination that enters in the integrand, 

will be called the microscopic spectral density^". This limit exists if chiral symmetry 
is broken. Our conjecture is that ps{u) is a universal function that only depends on 
the global symmetries of the QCD partition function. Because of universality it can 
be derived from the simplest theory with the global symmetries of the QCD partition 
function. Such theory is a chiral Random Matrix Theory which will be introduced later 
in these lectures. 

We emphasize again that the ?7a(1) symmetry of the QCD Dirac spectrum leads 
to two different types of eigenvalue correlations: spectral correlations in the bulk of 
the spectrum and spectral correlations near zero virtuality. The simplest example of 
correlations of the latter type is the microscopic spectral density defined in (|35|) . 

We close this subsection with two unrelated side remarks. First, the QCD Dirac 
operator is only determined up to a constant matrix. We can exploit this freedom 
to obtain a Dirac operator that is maximally symmetric. For example, the Wilson 
lattice QCD Dirac operator, D^, is neither Hermitean nor anti-Hermitean, but 
is Hermitean. 

Second, the QCD partition function can be expanded in powers of m? before or 
after averaging over the gauge field configurations. In the latter case one obtains sum 
rules for the inverse zeros of the partition function. As an example we quote. 



y- 

ml 



= (36) 



where we have averaged over field configurations with zero topological charge. 



SPECTRAL CORRELATIONS IN COMPLEX SYSTEMS 

Statistical Analysis of Spectra 

Spectra for a wide range of complex quantum systems have been studied both 
experimentally and numerically (a excellent recent review has been given by Guhr, 
Miiller-Groeling and Weidenmiiller^^). One basic observation is that the scale of varia- 
tions of the average spectral density and the scale of the spectral fluctuations separate. 
This allows us to unfold the spectrum, i.e. we rescale the spectrum in units of the local 
average level spacing. Specifically, the unfolded spectrum is given by 

= r{p{X'))dX', (37) 

J 00 

with unfolded spectral density 

Punf(A) = E^(^-ArO- (38) 

k 

The fluctuations of the unfolded spectrum can be measured by suitable statistics. 
We will consider the nearest neighbor spacing distribution, P{S), and moments of the 
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number of levels in an interval containing n levels on average. In particular, we will 
consider the number variance, S2(n), and the first two cumulants, 71 (n) and 72 (^^)- 
Another useful statistic is the A3(n)-statistic introduced by Dyson and Mehta"^^. It 
is related to T,2{n) via a smoothening kernel. The advantage of this statistic is that 
its fluctuations as a function of n are greatly reduced. Both E2(n) and A^ln) can be 
obtained from the pair correlation function defined as 

y2(A, A') = -(Punf (A)punf (AO) + (Punf (A)) (punf (A')) (39) 

Analytical expressions for the above statistics can be obtained for the eigenvalues 
of the invariant random matrix ensembles. They are defined as ensembles of Hermitean 
matrices with Gaussian independently distributed matrix elements, i.e. with probability 
distribution given by 

P{H) ~ e-"^^^'^. (40) 

Depending on the anti-unitary symmetry, the matrix elements are real, complex or 
quaternion real. They are called the Gaussian Orthogonal Ensemble (GOE), the Gaus- 
sian Unitary Ensemble (GUE) and the Gaussian Symplectic Ensemble (GSE), respec- 
tively. Each ensemble is characterized by its Dyson index (3 which is defined as the 
number of independent variables per matrix element. For the GOE, GUE and the GSE 
we thus have P — 1, 2 and 4, respectively. 

Independent of the value of /3, the average spectral density is the semicircle, 

N 



(p(A)) = -V2^. (41) 

71 

Analytical results for all spectral correlation functions have been derived for each 
of the three ensembles^^ via the orthogonal polynomial method. We only quote the 
most important results. The nearest neighbor spacing distribution, which is known 
exactly in terms of a power series, is well approximated by 

PiS)^ S'cxpi-aaS^), (42) 

where a^g is a constant of order one. The asymptotic behaviour of the pair correlation 
functions is given by^^ 

y2(A,A0- ^,^/_^,^, for = (43) 

1^2(A,A0~^g^^ for p = 2, (44) 

cos27r(A-A0 1 + (7r/2) sin27r(A - AQ 
^^^-^'^^ 4(A-A0 + An^iX - \'Y ^ " ^^^^ 

The 1/(A — A')^ tail of the pair correlation function results in a logarithmic dependence 
of the asymptotic behavior of ^2(71) and A3(n), 

E2(n) ~ (2/77^/3) logn and A3(n) ~ /3E2(n)/2. (46) 

Characteristic features of random matrix correlations are level repulsion at short dis- 
tances and a strong suppression of fluctuations at large distances. 

For uncorrelated eigenvalues the level repulsion is absent and one finds 

P(5)=exp(-5), (47) 

and 

S2(n) = n and A3(n) = n/15. (48) 
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Spectral Universality 

The main conclusion of numerous studies of eigenvalue spectra of complex systems 
is that spectral correlations of classically chaotic systems are given by RMT^^. As illus- 
tration we mention three examples from completely different fields. The first example is 
the nuclear data ensemble in which the above statistics are evaluated by superimposing 
level spectra of many different nuclei^^. The second example concerns correlations of 
acoustic resonances in irregular quartz blocks^^. In both cases the statistics that were 
considered are, within experimental accuracy, in complete agreement with the GOE 
statistics. The third example pertains to the zeros of Riemann's zeta function. Exten- 
sive numerical calculations^^ have shown that asymptotically, for large imaginary part, 
the correlations of the zeros are given by the GUE. 

The Gaussian random matrix ensembles introduced above can be obtaincd"^^ from 
two assumptions: i) The probability distribution is invariant under unitarity trans- 
formations; ii) The matrix elements are statistically independent. If the invariance 
assumption is dropped it can be shown with the theory of free random variables^^ 
that the average spectral density is still given by a semicircle if the variance of the 
probability distribution is finite. For example, if the matrix elements are distributed 
according to a rectangular distribution, the average spectral density is a semicircle. On 
the other hand, if the independence assumption is released the average spectral density 
is typically not a semicircle. For example, this is the case if the quadratic potential 
in the probability distribution is replaced by a more complicated polynomial potential 
V{H). Using the supersymmetric method for Random Matrix Theory, it was shown by 
Hackenbroich and Weidenmiiller^^ that the same supersymmetric nonlinear cr-model is 
obtained for a wide range of potentials V{H). This implies that spectral correlations 
of the unfolded eigenvalues are independent of the potential. Remarkably, this result 
could be proved for all three values of the Dyson index. 

An explicit construction of the correlation functions using orthogonal polynomials 
could only be performed^°' for (3 = 2. Several examples have been considered 

where both the invariance assumption and the independence assumption are relaxed. 
We mention H ^ H + A, where A is an arbitrary fixed matrix, and the probability 
distribution of H is given by a polynomial V{H). It was shown by P. Zinn-Justin^^ 
that also in this case the spectral correlations are given by the invariant random matrix 
ensembles. For a Gaussian probability distribution the proof was given by Brezin and 
Hikami^^. 

The domain of universality has been extended in the direction of real physical 
systems by means of the Gaussian embedded ensembles^^' The simplest example 
is the ensemble of matrix elements of n-particic Slater determinants of a two-body 
operator with random two-particle matrix elements. It can be shown analytically that 
the average spectral density is a Gaussian. However, according to substantial numerical 
evidence, the spectral correlations are in complete agreement with the invariant random 
matrix ensembles^^. 

Universal spectral correlations are obtained in the thermodynamic (or semi-classical) 
limit, oo, with {E — E')N fixed. Alternatively, one can take the thermody- 

namic limit with E — E' fixed. This leads to the Ambjorn-Jurkicwicz-Makeenko (AJM- 
) universality^^ for smoothened correlation functions^^. They are obtained from the ex- 
act correlation functions by replacing the oscillating terms by their average over a scale 
much larger than but small compared to the secular variation of the average spec- 
tral density. The result for the two-point correlation function is given by a smoothening 
of the leading order asymptotic result^^. However, correlations at macroscopic distance 
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in energy should include another ingredient. Namely, correlations resulting from the 
compactness of the support of the spectrum. Indeed, Ambjorn, Jurkiewicz and Ma- 
keenko showed^^ that the smoothened correlation functions for an arbitrary potential 
are determined completely by the endpoints of the spectrum. This theorem was proved 
for all three ensembles by Beenakker^° and it was generalized to multi-cut potentials by 
Akemann and Ambjorn^^. A general expression for the smoothened correlation func- 
tion in terms of one-point Greens functions has been derived for different deformations 
of the invariant random matrix ensembles^^' An interesting question is whether 

correlations given by this general expression are also found in real physical systems. 
There are several indications that the answer to this question is negative. First of all, 
AJM universality is closely related to the compactness of the support of the spectrum. 
In real physical systems, the spectral density usually increases with energy so that the 
average resolvent does not even exist (only the difference of two resolvents enters in 
the general expression). Second, in fully chaotic systems, with microscopic correlations 
given by Random Matrix Theory, it was shown by Berry^^ that the asymptotics of the 
two-point correlation function, as measured by the A3-statistic, is determined by the 
shortest periodic trajectory. 

Smoothened correlators are obtained via a perturbative expansion of the correla- 
tion functions. In Random Matrix Theory this is equivalent to a loop expansion in 1/N . 
The full non-perturbative result cannot be obtained this way. That requires the use of 
orthogonal polynomials or the super-symmetric method of Random Matrix Theory. 

Another type of universal behavior is given by correlations in the neighborhood of 
the largest eigenvalue. It was shown by Kanzieper and Preilikher^^ that for an arbitrary 
potential the spectral correlations at the soft edge of the spectrum are given by the Airy 
kernel. However, we are not aware of any physical applications of such correlations. 
Certainly, this type of universality is restricted to systems with a resolvent that is 
defined by a finite cuts on the real axis. For example, for the embedded ensembles with 
a Gaussian spectral density, we do not expect such type of universality. 

We have seen a large number of examples that fall into one of universality classes 
of the invariant random matrix ensembles. This calls out for a more general approach. 
Naturally, one thinks in terms of the renormalization group. This approach was pio- 
neered by Brezin and Zinn-Justin^^. The idea is to integrate out rows and columns of a 
random matrix and to show that the Gaussian ensembles are a stable fixed point. This 
was made more explicit in a paper by Higuchi et al.^. However, much more work is 
required to arrive at a natural proof of spectral universality. 

Although the above mentioned universality studies provide support for the validity 
of the Bohigas conjecture, the ultimate goal is to derive it directly from the underlying 
classical dynamics. An important first step in this direction was made by Berry^^. He 
showed that the asymptotics of the two-point correlation function is related to sum- 
rules for isolated classical trajectories. Another interesting approach was introduced 
by Andreev et al.^'^ who where able to obtain a supersymmetric nonlinear sigma model 
from spectral averaging. In this context we also mention the work of Altland and 
Zirnbauer^^ who showed that the kicked rotor can be mapped onto a supersymmetric 
sigma model. 
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CHIRAL RANDOM MATRIX THEORY 



Introduction of the Model 

In this section we will introduce an instanton liquid^^' ™ inspired chiral RMT for 
the QCD partition function. In the spirit of the invariant random matrix ensembles 
we construct a model for the Dirac operator with the global symmetries of the QCD 
partition function as input, but otherwise Gaussian random matrix elements. The 
chRMT that obeys these conditions is defined by^°' 



Nf 2 

Z^^,,(mi,.-.,m^,)= /w^det(P + m;)e-"^^^'^ (49) 

•' /=i 



where 



and W is a. n X m matrix with i/ — \n — m.\ and A'^ = n + m. The matrix elements of 
W are either real {P — 1, chiral Gaussian Orthogonal Ensemble (chGOE)), complex 
{f3 — 2, chiral Gaussian Unitary Ensemble (chGUE)), or quaternion real (/3 — 4. chiral 
Gaussian Symplectic Ensemble (chGSE)). As is the case in QCD, we assume that v 
does not exceed Vn, so that, to a good approximation, n = N/2. 

This model reproduces the following symmetries of the QCD partition function: 



The Ua{^) symmetry. All nonzero eigenvalues of the random matrix Dirac oper- 
ator occur in pairs ±A or are zero. 



• The topological structure of the QCD partition function. The Dirac matrix has 
exactly |z/| = |n — m| zero eigenvalues. This identifies i/ as the topological sector 
of the model. 



The flavor symmetry is the same as in QCD. For = 2 it is SU{Nf) x SU{Nf), 
for /? = 1 it is SU{2Nf) and for /? = 4 it is SU{Nf). 



• The chiral symmetry is broken spontaneously with chiral condensate given by 

E = lim 7ip(0)/N. (51) 

N—*oo 

{N is interpreted as the (dimensionless) volume of space time.) The symme- 
try breaking pattern is^=^ SU{Nf) x SU{Nf)/ SU{Nf), SU{2Nf) / Sp{Nf) and 
SU{Nf)/0{Nf) for /? = 2, 1 and 4, respectively the same as in QCD^^. 

• The anti-unitary symmetries. For three or more colors with fundamental fcrmions 
the Dirac operator has no anti-unitary symmetries, and gcncrically, the matrix 
elements of the Dirac operator arc complex. The matrix elements Wm of the 
corresponding random matrix ensemble are chosen arbitrary complex as well (/3 — 
2). For Nc — 2, the Dirac operator in the fundamental representation satisfies 

[CT2K,i'yD] = 0, (52) 

where C is the charge conjugation matrix and K is the complex conjugation 

operator. Because, {CT2KY = 1, the matrix elements of the Dirac operator can 
always be chosen real, and the corresponding random matrix ensemble is defined 
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with real matrix elements {j3 = 1). For two or more colors with gauge fields in the 
adjoint representation the anti-unitary symmetry of the Dirac operator is given 
by 

[CK,t-fD] = 0. (53) 

Because [CKY = —1, it is possible to rearrange the matrix elements of the Dirac 
operator into real quaternions. The matrix elements Wm of the corresponding 
random matrix ensemble are chosen quaternion real (/5 = 4). 

Together with the invariant random matrix ensembles, the chiral ensembles are 
part of a larger classification scheme. Apart from the random matrix ensembles dis- 
cussed in this review, this classification also includes random matrix models for disor- 
dered super-conductors^^. As pointed out by Zirnbauer^^, all known universality classes 
of Hermitean random matrices are tangent to the large classes of symmetric spaces 
in the classification given by Cartan. There is a one-to-one correspondence between 
this classification and the classification of the large families of Riemannian symmetric 
superspaces™ . 

Selected Results for the Chiral Random Matrix Ensembles 

The joint eigenvalue distribution follows from (^) by choosing the eigenvalues and 
eigenvectors as new integration variables. For Nj fiavors and topological charge u it is 
given by^^ 

P,(Ai, ■ ■ ■ , A„) = n - 1^ n Xr^^'-^'-' exp (-^ E ^l] ^ (54) 

k<l k \ k J 

where Cp^n is a normalization constant and z/ > 0. For (3 = 2 the average spectral 
density and the spectral correlation functions can be derived from (Q) with the help 
of the orthogonal polynomial method"^^. The associated polynomials are the general- 
ized Laguerre polynomials. That is why this ensemble is also known as the Laguerre 
ensemble^^' The spectral density and the two-point correlation function were also de- 
rived within the framework of the supersymmetric method of Random Matrix Theory™. 
The calculation of the average spectral density and the spectral correlations functions 
for (3=1 and /5 = 4 is much more complicated. However, with the help of skew- 
orthogonal polynomials®^' ®^ exact analytical results for finite N can be obtained as 
well. 

From the properties of the Laguerre polynomials it can be shown that, independent 
of the value of (3, the average spectral density is a semi-circle 



p(A) = (nSV7r)y'4/S2 - A^. (55) 



The microscopic spectral density can be derived from the limit ([35| ) of the exact 
spectral density for finite A^. For = 3, Nf fiavors and topological charge v it is given 



PS{Z) = ^ [JI{Z) - Ja+l{z)Ja-l{z)) , (56) 

where a = Nj + |z^|. The expressions for SU{2) with fundamental fermions {(3 = 1) are 
much more complicated. In this case we find the microscopic spectral density®^ 



Ps{z) = —J2a+l{zT) + — / dw{zw) 

4 z JO 



w dw z dz 



wJ2a{zE)J2a-l{wE) - Z J2a-l{zE) J2a{wT,) 
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where a is the combination 

a = Nf-^- + ^. (58) 

The microscopic spectral density in the symmetry class with (3 = 4 was first calculated 
by Nagao and Forrester It is given by 

Ps{z) = 2z^ duU^ dv [JAa-l{'^UVZ)J4,a{2uz)—vJAa-l{'^UZ)J4,a{2uVz)\ (59) 

JO Jo 

with Aa = Nf + 2\v\ + 1. 

The spectral correlations in the bulk of the spectrum are given by the invariant 
random matrix ensemble with the same value of (3. For (3 = 2 this was already shown 
three decades ago by Fox and Kahn'''^. For (3 = 1 and /5 = 4 this was only proved 
recently^^. 

Duality between Flavor and Topology 

As one can observe from the joint eigenvalue distribution, for (3 = 2 the dependence 
on Nf and v enters only through the combination Nf + u. This allows for the possibility 
of trading topology for fiavors. In this section we will work out this duality for the 
finite volume partition function. This relation completes the proof of the conjectured 
expression^^ for the finite volume partition function for different quark masses and 
topological charge z/. 

In terms of the eigenvalues the partition function ( ^91 ) is given by 



k f f,k 

(60) 

where the Vandermonde determinant is defined by 

A{Xl)=Yl{Xl-X^). (61) 

k,l 

By inspection we have 

^iVy,v(mi,---,mjv^) ^ . n feio\ 

ZNf+ufiimi,- ■ ■ ,mNf,0,- ■ ■ ,0), (62) 



where the argument of the last factor has u zeros. 

As an example, the simplest nontrivial identity of this type is given by 

Zi^i(m) ~ mZ2fl{m, 0). (63) 

Let us prove this identity without relying on Random Matrix Theory. According to the 
definition (pHj) we have 

Zi,i{m) ~ j dee''e"'^^''°'\ (64) 

and 

Ju&(2) 
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where M is a diagonal matrix with diagonal elements m and 0. The integral over U 
can be performed by diagonalizing U according to U = Uie^'^'^Ui^ , and choosing Ui 
and (pk as new integration variables. The Jacobian of this transformation is 

J ~ A^e"'"'). (66) 

The integral over Ui can be performed using the Itzykson-Zuber formula. This results 



m 



J m(cos0i — cos<p2) 

Both terms in the last factor result in the same contribution to the integral. Let us 
consider only the first term ~ exp(ml^S cos0i). Then the integral over 02 has to be 
defined as a principal value integral. If we use the identity 

_g-ife|2 / ^sin0^^\ 

7 ; — = / COS 01 — , , , , (68) 

(cos(/.i -COS02) V tan((0i + 02)/2)y ^ ' 

the 02-integral of the term proportional to sin 01 gives zero because of the principal 
value prescription. The term proportional to cos0i trivially results in Zi i. We leave 
it as an exercise to the reader to generalize this proof to arbitrary Nf and u. 

The group integrals in finite volume partition function (^SD were evaluated by 
Leutwyler and Smilga* for equal quark masses. An expression for different quark 
masses was obtained by Seneret al.^^. The expression could only be proved for u = 0. 
The above duality can be used to relate a partition function at arbitrary u to a partition 
function at = 0. This completes the proof of the conjectured expression for arbitrary 
topological charge. 



UNIVERSALITY IN CHIRAL RANDOM MATRIX THEORY 

In the chiral ensembles, two types of universality studies can be performed. First, 
the universality of correlations in the bulk of the spectrum. As discussed above, they 
are given by the invariant random matrix ensembles. Second, the universality of the 
microscopic spectral density and the eigenvalue correlations near zero virtuality. The 
aim of such studies is to show that microscopic correlations are stable against deforma- 
tions of the chiral ensemble away from the Gaussian probability distribution. Recently, 
a number of universality studies on microscopic correlations have appeared. They will 
be reviewed in this section. 

In a first class of universahty studies one considers probability distributions that 
maintain unitary invariance, i.e. 

y2 oo 

P{W) ~ exp{-N(3— akTriW^W)''). (69) 
^ k=i 

The first study of this kind was performed by Brezin, Hikami and Zee^^. They consid- 
ered a potential with only oi and a2 different from zero and showed that the microscopic 
spectral density is independent of 02. A general proof valid for arbitrary potential was 
given by Nishigaki^^' The essence of the proof is a remarkable generalization of the 
identity for the Laguerre polynomials, 

lim L„(-) = Jo(2v^) , (70) 
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to orthogonal polynomials determined by an arbitrary potential V. This relation was 
proved by deriving a differential equation from the continuum limit of the recursion 
relation for orthogonal polynomials. In a consecutive work, Akemann, Damgaard, 
Magnea and Nishigaki^^ extended this proof to all microscopic correlation functions. 

In a second class of universality studies one considers deformations of the Gaus- 
sian random matrix ensemble that violate unitary invariance. In particular, one has 
considered the case where the matrix W in (^) is replaced by 

W-^W + A, (71) 

where, because of the unitary invariance, the matrix A can always be chosen diagonal. 
The simplest case with A = ttT times the identity was considered by Sener et al}^. This 
model provides a schematic model of the chiral phase transition. For large matrices, 
the average resolvent defined by 

G{z) = Tr( ^—-) (72) 



obeys the cubic equation®^' (the parameter S = 1 in (^)) 

- 2zG^ + G{z^ - TT^T^ + l)-z = 0. (73) 

This equation was first obtained using a diagrammatic method®^' Later, this deriva- 
tion was rewritten*^ in terms of the so called blues function^^. The average spectral 
density given by 

p(A) = --lmG{z = A) (74) 

TT 

is a semicircle at vrT = and splits into two arcs at ttT = 1. For the spectral density at 
zero one obtains p(0) = \/l — vr^T^, and therefore chiral symmetry is broken for irT < 1, 
and is restored above this temperature. In spite of this drastic change in average spectral 
density, it could be shown^^ with the help of a supersymmetric formulation of Random 
Matrix Theory that the microscopic spectral density does not depend on T. 

The super-symmetric method in the first paper by Sener et al}^ is not easily 
generalizable to higher order correlation functions. A natural way to proceed is to 
employ the super-symmetric method introduced by Guhr^^. In the case of /3 = 2 this 
method results in an expression for the kernel determining all correlation functions. 
This approach was followed in two papers, one by Guhr and Wettig^^ and one by 
Sener et al}^. The latter authors studied microscopic correlation functions for A in 
( |7ID proportional to the identity, whereas Guhr and Wettig considered an arbitrary 
diagonal matrix A. It was shown that independent of the matrix A, the correlations 
are given by the Bessel kernel^^. Of course, a necessary condition on the matrix A is 
that chiral symmetry is broken. Guhr and Wettig also showed that correlations in the 
bulk of the spectrum are insensitive to A. 

The main ingredient of the proof was a supersymmetric generalization^^ of the 
Itzykson-Zuber-Harish-Chandra type integral^^' 

I dfi{U, V) exp(ReTrt/t^l^i?) = C ^ ^^^kAUnsi)) 
^ ^ ^ ^ ^ nfli(rfcS,)'^A(52)A(i?2) ^ ^ 

Here, U e U{Ni), V G f/(A^2), = Ni - N2 > 0. The positive square roots of the 
eigenvalues of S'^S and i?-R^ are denoted by, Sk and r^, respectively. The integral is 
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over the invariant measure. The constant C in the r.h.s. can be evaluated, and tlie 
are modified Bessel functions. 

The deformation W —>■ W + A with the probabihty distribution for W given by 
an arbitrary invariant potential has not yet been considered. We have no doubt that 
universality proofs along the lines of methods developed by P. Zinn-Justin^^ can be 
given. 

We wish to emphasize that all universality studies for the chiral ensembles have 
been performed for j3 = 2. The reason is that /? = 1 and /? = 4 are mathematically 
much more involved. It certainly would be of interest to extend the above results to 
these cases as well. 

In addition to the above analytical studies the universality of the microscopic 
spectral density also follows from numerical studies of models with the symmetries of the 
QCD partition function. In particular, we mention strong support in favor universality 
from a different branch of physics, namely from the theory of universal conductance 
fluctuations. In that context, the microscopic spectral density of the eigenvalues of the 
transmission matrix was calculated for the Hofstadter^^ model, and, to a high degree 
of accuracy, it agrees with the random matrix prediction^^. Other studies deal with a 
class random matrix models with matrix elements with a diverging variance. Also in 
this case the microscopic spectral density is given by the universal expressions^'' (|57| ) 
and (H). 

The conclusion that emerges from all numerical and analytical work on modifled 
chiral random matrix models is that the microscopic spectral density and the correla- 
tions near zero virtuality exhibit a strong universality that is comparable to the stability 
of microscopic correlations in the bulk of the spectrum. 

Of course, QCD is much richer than chiral Random Matrix Theory. One question 
that should be asked is at what scale (in virtuality) QCD spectral correlations deviate 
from RMT. This question has been studied by means of instanton liquid simulations. 
Indeed, at macroscopic scales, it was found that the number variance shows a linear 
dependence instead of the logarithmic dependence observed at microscopic scales^^. 
More work is needed to determine the point where the crossover between these two 
regimes takes place. 



LATTICE QCD RESULTS 

Recently, the Dirac spectrum in lattice QCD received a good deal of attention. In 
particular, the connection between the topology of fleld conflgurations and the spectrum 
of the Wilson Dirac operator has been studied in detail^^' Other studies are 

related to the connection between the Wilson Dirac spectrum and the localization 
properties of the eigenfunctions^^^ . 

In this section we will focus ourselves on the spectral correlations of the lattice QCD 
Dirac operator. Both correlations in the bulk of the spectrum and the microscopic 
spectral density will be studied. Consistent with universality arguments presented 
above, we flnd that spectral correlations are in complete agreement with chiral Random 
Matrix Theory. 

Correlations in the Bulk of the Spectrum 

Recently, Kalkreuter^^ calculated all eigenvalues of the Nc = 2 lattice Dirac oper- 
ator both for Kogut-Susskind (KS) fermions and Wilson fermions for lattices as large 
as 12^. For the Kogut-Susskind Dirac operator, D^^ , we use the convention that it is 
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anti-Hermitean. Because of the Wilson-term, the Wilson Dirac operator, D , is neither 



Hermitean nor anti-Hermitean. Its Hermiticity relation is given by D 



Therefore, the operator 75/^^ is Hermitean. However, it does not anti-commute with 
75, and its eigenvalues do not occur in pairs ±Afe. 
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Figure 2. Spectral correlations of Dirac eigenvalues for Wilson fermions (left) and KS-fermions 
(right). Results are shown for the number variance, ^2{n), the A3— statistic and the nearest 
neighbor spacing distribution, P{S). The full, dashed and dotted curves represent the analytical 
result for the GOE, GUE and GSE, respectively. 



In the the case of SU{2), the anti-unitary symmetry of the Kogut-Susskind and the 
Wilson Dirac operator is given by^°^' 

[L'^^,T2X] = 0, and [75L'^,75CXt2] =0. (76) 

Because 

{T2Kf^-l, and (75CXt2)' = 1, (77) 

the matrix elements of the KS Dirac operator can be arranged into real quaternions, 
whereas the Wilson Dirac operator can be expressed into real matrix elements. There- 
fore, we expect that eigenvalue correlations in the bulk of the spectrum are described by 
the GSE and the GOE, respectively^^. The microscopic correlations for KS fermions are 
described by the chGSE. However, the microscopic correlations for Wilson fermions are 
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not described by the chGOE but rather by the GOE. Because of the anti-unitary sym- 
metry, the eigenvalues of the KS Dirac operator are subject to the Kramers degeneracy, 
i.e. they are double degenerate. 

In both cases, the Dirac matrix is tri-diagonalized by CuUum's and Willoughby's 
Lanczos procedure^^^ and diagonalized with a standard QL algorithm. This improved 
algorithm makes it possible to obtain all eigenvalues. This allows us to test the accuracy 
of the eigenvalues by means of sum-rules for the sum of the squares of the eigenvalues 
of the lattice Dirac operator. Typically, the numerical error in the sum rule is of order 
10-8. 
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Figure 3. The number variance, T,2{n) and tiie first two cumulants, 71 (n) and 72(^1) as a function of 
n for eigenvalues of the Wilson Dirac operator (left) and the Kogut-Susskind Dirac operator (right). 
The full and dotted curves represent the analytical result for the GOE and the GSE, respectively. 



As an example, in Fig. 1 we show a histogram of the overall Dirac spectrum for 
KS fermions a.t (3 = 2.4. Results for the spectral correlations are shown in Figs. 2 and 
3. The results for KS fermions are for 4 dynamical flavors with ma — 0.05 on a 12^ 
lattice. The results for Wilson fermions were obtained for two dynamical flavors on 
a 8^ X 12 lattice. For the values of j3 and k we refer to the labels of the figure. For 
13 > 2.4, with our lattice parameters for KS fermions, the Dirac spectrum near zero 
virtuality develops a gap. Of course, this is an expected feature of the weak coupling 
domain. For small enough value of k, the Wilson Dirac spectrum shows a gap at A = 
as well. In the scaling domain the value of k is just above the critical value of k. A 
qualitative description of the Wilson Dirac spectrum can be obtained with a random 
matrix model with the structure of the Wilson Dirac operator^°^ . 
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The eigenvalue spectrum is unfolded by fitting a second order polynomial to the 
integrated spectral density of a stretch of 500-1000 eigenvalues. The results for S2(n), 
A^^n) and P{S) in Fig. 2 show an impressive agreement with the RMT predictions. 
The fluctuations in ^2(77.) are as expected from RMT. The advantage of As-statistic is 
well illustrated by this figure. We also investigated^^ the n dependence of the first two 
cumulants of the number of levels in a stretch of length n. Results presented in Fig. 4 
show a perfect agreement with RMT. 

Spectra for different values of /3 have been analyzed as well. It is probably no 
surprise that random matrix correlations are found at stronger couplings. What is 
surprising, however, is that even in the weak-coupling domain (/3 = 2.8) the eigenvalue 
correlations are in complete agreement with Random Matrix Theory. Finally, we have 
studied the stationarity of the ensemble by analyzing level sequences of about 200 eigen- 
values (with relatively low statistics). No deviations from random matrix correlations 
were observed all over the spectrum, including the region near A = 0. This justifies the 
spectral averaging which results in the good statistics in Figs. 2 and 3. 

In the case of three or more colors with fundamental fermions, both the Wilson and 
Kogut-Susskind Dirac operator do not possess any anti-unitary symmetries. Therefore, 
our conjecture is that in this case the spectral correlations in the bulk of the spectrum 
of both types of fermions can be described by the GUE. In the case of two fundamental 
colors the continuum theory and Wilson fermions are in the same universality class. 
It is an interesting question of how spectral correlations of KS fermions evolve in the 
approach to the continuum limit. Certainly, the Kramers degeneracy of the eigenvalues 
remains. However, since Kogut-Susskind fermions represent 4 degenerate fiavors in the 
continuum limit, the Dirac eigenvalues should obtain an additional two-fold degeneracy. 
We are looking forward to more work in this direction. 

The Microscopic Spectral Density 

The advantage of studying spectral correlations in the bulk of the spectrum is 
that one can perform spectral averages instead of ensemble averages requiring only a 
relatively small number of equilibrated configurations. This so called spectral ergodicity 
cannot be exploited in the study of the microscopic spectral density. In order to gather 
sufficient statistics for the microscopic spectral density of the lattice Dirac operator a 
large number of independent configurations is needed. One way to proceed is to generate 
instanton-liquid configurations which can be obtained much more cheaply than lattice 
QCD configurations. Results of such analysis^"^ show that for Nc = 2 with fundamental 
fermions the microscopic spectral density is given by the chGOE. For A^c = 3 it is given 
by the chGUE. One could argue that instanton-liquid configurations can be viewed 
as smoothened lattice QCD configurations. Roughening such configurations will only 
improve the agreement with Random Matrix Theory. 

Of course, the ultimate goal is to test the conjecture of microscopic universality 
for realistic lattice QCD configurations. In order to obtain a very large number of 
independent gauge field configurations one is necessarily restricted to relatively small 
lattices. The first study in this direction was reported recently^^' In this work, the 
quenched SU{2) Kogut-Susskind Dirac operator is diagonalized for lattices with linear 
dimension of 4, 6, 8 and 10, and a total number of configurations of 9978, 9953, 3896 
and 1416, respectively. The results were compared with predictions from the chGSE. 

We only show results for the largest lattice. For more detailed results, including results 
for the two-point correlation function, we refer to the original work. In Fig. 4 we show 
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Figure 4. The distribution of the smallest eigenvalue (left) and the microscopic spectral density 
(right) of the Kogut-Susskind Dirac operator for two colors and (3 — 2.0. Lattice results are 
represented by the histogram, and the analytical results for the chGSE are given by the dashed 
curves. 



the distribution of the smallest eigenvalue (left) and the microscopic spectral density 
(right). The lattice results are given by the full line. The dashed curve represents 
the random matrix results. The distribution of the smallest eigenvalue was derived by 
Forrester^"^ and is given by 

P(A„.in) = «y|(«A^i,)3/2j3/2(«A^i,)e-^(°^-'")', (78) 

where a = VTj. The random matrix result for the microscopic spectral density is given 
in eq. (|59|). We emphasize that the theoretical curves have been obtained without any 
fitting of parameters. The input parameter, the chiral condensate, is derived from the 
same lattice calculations. The above simulations were performed at a relatively strong 
coupling of /3 = 2. Recently, the same analysis^^^ was performed for f3 = 2.2 and for 
f3 = 2.5 on a 16^ lattice. In both cases agreement with the random matrix predictions 
was found^°^. 

An alternative way to probe the Dirac spectrum is via the valence quark mass 
dependence of the chiral condensate^^ defined as 

E(™) = l/<(A(p(A))^. (79) 

The average spectral density is obtained for a fixed sea quark mass. For masses well 
beyond the smallest eigenvalue, S(m) shows a plateau approaching the value of the 
chiral condensate S. In the mesoscopic range (|l]), we can introduce u = XmN and 
X = mNTj as new variables. Then the microscopic spectral density enters in E(m). For 
three fundamental colors the microscopic spectral density for (3 = 2 (eq. (|56|)) applies 
and the integral over A in ( |75D can be performed analytically. The result is given by^^, 

^ = x{Ia{x)Ka{x) + Ia+l{x)Ka-l{x)), (80) 

where a = Nf + In Fig. 2 we plot this ratio as a function of x (the 'volume' V is 
equal to the total number of Dirac eigenvalues) for lattice data of two dynamical flavors 
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with mass ma = 0.01 and A^c = 3 on a 16^ x 4 lattice. We observe that the lattice data 
for different values of (3 fall on a single curve. Moreover, in the mesoscopic range this 
curve coincides with the random matrix prediction for Nf = u = 0. 
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Figure 5. The valence quark mass dependence of the chiral condensate S^(m) plotted as S]^(m)/E 
versus mVT,. The dots and squares represent lattice results by the Columbia group^^ for values of (3 
as indicated in the label of the figure. 



Apparently, the zero modes are completely mixed with the much larger number of 
nonzero modes. For eigenvalues much smaller than the sea quark mass, one expects 
quenched {Nf = 0) eigenvalue correlations. In the same figure the dashed curves 
represent results for the quark mass dependence of the chiral condensate (i.e. the mass 
dependence for equal valence and sea quark masses). In the sector of zero topological 
charge one finds^^°' ^' ® 



S^) hiu) 
E 

and 



S Io{u) 



u{PM-ll{u)) 



for Nf = 1, (81) 



for Nf = 2. (82) 



We observe that both expressions do not fit the data. Also notice that, according 
to Gockeler et al}^^, eq. ( ^T]) describes the valence mass dependence of the chiral 
condensate for non-compact QED with quenched Kogut-Susskind fermions. However, 
we were not able to derive their result (no derivation is given in the paper). 



CHIRAL RANDOM MATRIX THEORY AT fi ^ 

At nonzero temperature T and chemical potential fi a schematic random matrix 
matrix model of the QCD partition function is obtained by replacing the Dirac operator 

in (H) by87' 112, 88, 35 

-n ( tW + tnT + fi\ , . 
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Here, Qt = T ®n (2ri + l)7rl + fi are the matrix elements of 70^0 + /i7o in a plane 
wave basis with anti-periodic boundary conditions in the time direction. Below, we 
will discuss a model with Qt absorbed in the random matrix and /i 7^ 0. The aim of 
this model is to explore the effects of the non-Hermiticity of the Dirac operator. For 
example, the random matrix partition function (EDf) with the Dirac matrix (RHf) is well 



suited for the study of zeros of this partition function in the complex mass plane and 
in the complex chemical potential plane. For a complete analytical understanding of 
the location of such zeros we refer to the work by Halasz et al}^^. 

The term /i7o does not affect the anti-unitary symmetries of the Dirac operator. 
This is also the case in lattice QCD where the color matrices in the forward time 
direction are replaced by f/ — > e^f/ and in the backward time direction by f/^ e~^f/^ 
For this reason the universality classes are the same as at zero chemical potential. 

The Dirac operator that will be discussed in this section is thus given by 

^W-f-W^/'o'^'V (84) 



where the matrix elements W are either real (/3 = 1), complex [jS = 2) or quaternion 
real (/3 = 4). For all three values of (3 the eigenvalues of are scattered in the 

complex plane. 

Since many standard random matrix methods rely on convergence properties based 
on the Hermiticity of the random matrix, direct application of most methods is not 
possible. The simplest way out is the Hermitization^^'^ of the problem, i.e we consider 
the Hermitean operator 

For example, the generating function in the supersymmetric method of Random Matrix 
Theoryi^S' is then given by^^^' 

^^J^J^ = ^ det(D^(.,.-)) 

The determinants can be rewritten as fermionic and bosonic integrals. Convergence 
is assured by the Hermiticity and the infinitesimal increment k. The resolvent follows 
from the generating function by differentiation with respect to the source terms 



(87) 

J=J*=0 



Notice that, after averaging over the random matrix, the partition function depends in 
a non-trivial way on both z and z* . The spectral density is then given by 

p{z) = --^G{z). (88) 
vr oz* 

Alternatively, one can use the replica trick^^^' with generating function given by 

{dei^fDH{z,z*)). (89) 

The idea is to perform the calculation for integer values of Nf and perform the limit 
Nf — i> at the end of the calculation. Although the replica limit, Nf — > 0, fails in 
general^^^, it is expected to work for det Dh{z, z*) because it is positive definite (or 
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zero). Then the partition function is a smooth function of Nf. For other techniques 
addressing nonhermitean matrices we refer to the recent papers by Feinberg and Zee^^^ 
and Nowak and co-workers^^^. One recent method that does not rely on the Hermitic- 
ity of the random matrices is the method of complex orthogonal polynomials'^^. This 
method was used by Fyodorov et al.^^ to calculate the number variance and the near- 
est neighbor spacing distribution in the regime of weakly nonhermitean matrices. As 
surprising new result, they found an S^^"^ repulsion law. 

In the physically relevant case of QCD with three colors, the fermion determinant 
is complex for nonzero chemical potential. Its phase prevents the convergence of fully 
unquenched Monte-Carlo simulations (see Kogut et al}'^'^ for the latest progress in this 
direction). However, it is possible to perform quenched simulations. In such calculations 
it was found that the critical chemical potential Hc ~ \/rn, instead of a third of the 
nucleon mass^^^. This phenomenon was explained analytically by Stephanov^^ with 
the help of the above random matrix model. He could show that for small /i the 
eigenvalues are distributed along the imaginary axis in a band of width ~ /i^ leading 
to a critical chemical potential of /ic ~ m^. As has been argued above, the quenched 
limit is necessarily obtained from a partition function in which the fermion determinant 
appears as 

lim |detP(/i)|^^ (90) 

instead of the same expression without the absolute value signs. The partition function 
with the absolute value of the determinant can be interpreted as a partition function of 
an equal number of fermions and conjugate fermions. The critical value of the chemical 
potential, equal to half the pion mass, is due to Goldstone bosons with net a baryon 
number consisting of a quark and conjugate quark. The reason that the quenched limit 
does not correspond to the standard QCD partition function is closely related to the 
failure of the replica trick in the case of a determinant with a nontrivial phase. 

Both for (3 = 1 and /3 = 4 the fermion determinant, det(T'(/i) + m), is real. This 
is obvious for (3=1. For (3 = 4: the reality follows from the identity q* = 0^(102 for a 
quaternion real element g, and the invariance of a determinant under transposition. We 
thus conclude that quenching works for an even number of flavors. Consequently, chiral 
symmetry will be restored for arbitrarily small nonzero yU, whereas a condensate of a 
quark and a conjugate quark develops. Indeed, this phenomenon has been observed in 
the strong coupling limit of lattice QCD with two colors^^^. 

In the quenched approximation, the spectral properties of the random matrix 
ensemble (0) can be easily studied numerically by simply diagonalizing a set of matrices 
with probability distribution (|^). In Fig. 6 we show numerical results^^^ for the 
eigenvalues of a few 100 x 100 matrices for = 0.15 and \i = 0.5. The dots represent 
the eigenvalues in the complex plane. The full line is the analytical result^^ for the 
boundary of the eigenvalues which is given by the algebraic curve 

This result was obtained by Stephanov for (3 = 2 with the help of the fermionic replica 
trick. This amounts to rewriting the determinants in ( PP] ) as integrals over Grassmann 
variables. Since Grassmann integrals are always convergent the infinitesimal increment 
K can be put equal to zero. This method can be extended^^^ to (3 = 1 and (3 = 4. 
Although the effective partition function is much more complicated, it can be shown 
without too much effort that the solutions of the saddle point equations are the same if 
the variance of the probability distribution is scaled as 1//3. In particular, the boundary 
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Figure 6. Scatter plot of the real (x), and the imaginary parts ( y) of the eigenvalues of the 
random matrix Dirac operator. The values of f3 and /i are given in the labels of the figure. The full 
curve shows the analytical result for the boundary. 



of the domain of eigenvalues is the same in each of the three cases. However, as one 
observes from Fig. 6, for /3 = 1 and (3 = 4 the spectral density deviates significantly 
from the saddle-point result. For /3 = 1 we find an accumulation of eigenvalues on the 
imaginary axis, whereas for /3 = 4 we find a depletion of eigenvalues in this domain. 
This depletion can be understood as follows. For = all eigenvalues are doubly 
degenerate. This degeneracy is broken at /i 7^ which produces the observed repulsion 
of the eigenvalues. 

The number of purely imaginary eigenvalues for (3 = 1 appears to scale as \/N. 
This explains that this effect is not visible in a leading order saddle point analysis. 
From a perturbative analysis of (|89|) one obtains a power series in 1/N. Clearly, the 
^ dependence requires a truly nonperturbative analysis of the partition function (|i9| ) 



with the Dirac operator (84). Such a. yN scaling behavior is typical for the regime of 
weak non-hermiticity first identified by Fyodorov et al}^^. Using the supersymmetric 



method for the generating function (^) the v dependence was obtained analytically 
by Efetov^20_ 

A similar cut below a cloud of eigenvalues was found in instanton liquid simulations^^® 
for Nc = 2 aX fi and in a random matrix model of arbitrary real matrices^^^. The 
depletion of the eigenvalues along the imaginary axis was observed earlier in lattice 
QCD simulations with staggered fermions^^". Obviously, more work has to be done in 
order to arrive at a complete characterization of the universal features^^ in the spectrum 
of nonhermitean matrices. 



CONCLUSIONS 

We have argued that there is an intimate relation between correlations of Dirac 
eigenvalues and the breaking of chiral symmetry. In the chiral limit, the fermion deter- 
minant suppresses gauge field configurations with small Dirac eigenvalues. Correlations 
counteract this suppression, and are a necessary ingredient of chiral symmetry break- 
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ing. From the study of eigenvalue correlations in strongly interacting systems, we have 
concluded that they are described naturally with by Random Matrix Theory with the 
global symmetries of the physical system. In QCD, this led to the introduction of 
chiral Random Matrix Theories. They provided us with an analytical understanding 
of the statistical properties of the eigenvalues on the scale of a typical level spacing. 
In particular, impressive agreement between lattice QCD and chiral Random Matrix 
Theory was found for the microscopic spectral density and for spectral correlations in 
the bulk of the spectrum. An extension of this model to nonzero chemical potential 
explains some intriguing properties of previously obtained lattice QCD Dirac spectra 
and instanton liquid Dirac spectra. 
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